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It will be shown that a given realization of nonlinear electrodynamics, used as source of Einstein's 
equations, generates a cosmological model with interesting features, namely a phase of current 
cosmic acceleration, and the absence of an initial singularity, thus pointing to a way to solve two 
important problems in cosmology. 



I. INTRODUCTION 

In spite of the success of the Standard Cosmologi- 
cal Model (SCM), several problems remain to be solved, 
most notably the cause of the current acceleration stage, 
and the initial singularity. The observation that the uni- 
verse is undergoing a phase of accelerated expansion has 
sparkled an intense activity directed to investigate the 
possible candidates to fuel the acceleration. One possi- 
bility is the modification of the gravitational action by 
the addition to the Einstein-Hilbert action of terms that 
depend on some power of the curvature. The simplest of 
these models |l| with action given by 



S 




d x. 



has been shown to disagree with solar system observa- 
tions Q. Later it was proved that more general models 
Hi [J i such as those given by the Lagrangian 



R 7 ' 



(1) 



can describe a phase of accelerated expansion (controlled 
by the negative powers of n) and produce other modifi- 
cations in the strong regime (where terms with n > are 
dominant). 

A second possibility relates the acceleration to the mat- 
ter sector of the theory, either by incorporating a cosmo- 
logical constant or by postulating the existence of mat- 
ter fields with peculiar properties. In the latter case, the 
most popular choice is a scalar field in the presence of a 
potential Q , or with a nontrivial Lagrangian (see [|| and 
references therein). Yet another possibility is a vector 
field, as discussed in Q. Although the addition of these 
fields to the action yields acceleration (and in some cases 
dark matter as a byproduct the need of (yet to be 
observed) matter with unusual properties is certainly a 
hindrance. 

Another problem in the description furnished by the 
SCM is that of the initial singularity. It is a well-known 



fact that under certain assumptions, the SCM unavoid- 
ably leads to a singular behaviour of the curvature in- 
variants in what has been termed the Big Bang. This is 
a highly distressing state of affairs, because the presence 
of a singularity entails the breakdown of all known phys- 
ical theories. There are hopes that quantum gravity will 
solve this deficiency but, since this theory is still incom- 
plete, one is obliged to explore alternative routes in order 
to avoid the initial singularity. 

We shall show in this article that nonlinear electrody- 
namics (NLED) can be useful in the discussion of possible 
solutions to these two problems of the SCM. Regarding 
the accelerated expansion, an alternative to the two pos- 
sibilities presented above was given in the model intro- 
duced in 0, where the action for the electromagnetic 
field was modified by the addition of a new term, namely 



(2) 



with F = F M „ . Through a suitable averaging process, 
this action yields an accelerated expansion phase for the 
evolution of the universe in the weak-EM-field regime. 
It also correctly describes the electric field of an isolated 
charge Q without the need of averaging. 

Scenarios that avoid the initial singularity have been 
intensely studied over the years. As examples of the lat- 
est realizations we can quote the pre-big-bang universe 
[Tfj| and the ekpyrotic universe 0. While these mod- 
els are based on deep changes to known physics (involv- 
ing new entities as scalar fields or branes) the model 
we present here relies instead on an unique entity, the 
electromagnetic field. The new ingredient that we in- 
troduce is a modification of the dynamics, which differs 
from that of Maxwell in certain regimes. Specifically, the 
Lagrangian we will work with is given by 




aF 2 



1 



F 



1 
F' 



(3) 



The dimensional constants a and 7 are to be deter- 
mined by observation. We shall see that in Friedmann- 
Lemaitre-Robertson- Walker (FLRW) geometry the first 



2 



term dominates in very early epochs, the Maxwell term 
dominates in the radiation era, and the last term is re- 
sponsible for the accelerated phase. 

It will be shown in this article that the Lagrangian 
given in Eqn. © yields a unified scenario to describe both 
the acceleration of the universe (for weak fields, as dis- 
cussed above) and the avoidance of the initial singular- 
ity as a consequence of its properties in the strong-field 
regime (see also ,$]). The plan of the article is as follows. 
In section II we present some consequences of applying 
Tolman average procedure to a nonlinear EM field con- 
figuration in a FLRW geometry, we describe its energy 
content and show how nonlinearities change the equa- 
tion of state of the electromagnetic field. In section III 
we present the notion of magnetic universe, analyze the 
related energy density and pressure and show the unex- 
pected result that the dependence of the magnetic field 
on the scale factor a(t) is the same irrespectively of the 
form of the Lagrangian. Section IV deals with the condi- 
tions of the existence of a bouncing and the acceleration 
for a generic fluid. In Section V we present a simple 
model of a bouncing universe dominated by the nonlin- 
ear properties of the averaged magnetic field, as well as 
a discussion about a cyclic universe and its correspond- 
ing properties in this model. Section VI presents the 
modifications that nonlinearities of the electromagnetic 
field in the weak regime produce for a late phase of the 
universe. Several interesting topics related to nonlinear 
electromagnetism (such as the field of a point charge and 
causality) are discussed in the appendices. 

Section VII presents a complete scenario in which both 
problems - bouncing and acceleration - are treated in a 
single and an unified theory. 



Following a standard procedure we define the vol- 
umetric spatial average of a quantity X at the time t 
by 



X 



-g d x, 



(6) 



where V = j yj—g d 3 x and Vq is a sufficiently large time- 
dependent three-volume. In this notation, for the elec- 
tromagnetic field to act as a source for the FLRW model 
we need to impose that 



Ei — 0, Hi — 0, EiHj — 0, 



1 , 1 , 



(7) 



EiEj = --E gtj, HiHj — ~-^H g tj . (8) 

With these conditions, the energy-momentum tensor of 
the EM field associated to C — C(F, G) can be written 
as that of a perfect fluid, 



Tftu = (fi + p)v^v v -pg^ v , 



(9) 



where 



P = -C + GC G -4C F E 2 , (10) 
p = t-G L G -\(2H 2 - E 2 )C F , (11) 



and Cj 



dC 

dA ■ 



A = F,G. In the present work we shall 



restrict to the case of a nonlinear theories defined by £ = 
C(F). In such a case, the energy-momentum is diagonal, 
with 



II. THE AVERAGE PROCEDURE AND THE 
FLUID REPRESENTATION 



The effects of a nonlinear electromagnetic theory in a 
cosmological setting have been studied in several articles 
|12|. In particular, it has been shown that under certain 
assumptions (as for instance the existence of a compact- 
ified space) a Born-Infeld field as a source can model 
the accelerated expansion of the universe [l^. although 
this is not possible in the original formulation of the the- 
ory, as we will show in Appendix. Notice that due to the 
isotropy of the spatial hyper-surfaces of FLRW geometry, 
an average procedure is needed if electromagnetic fields 
are to act as a source of gravity. Given a generic gauge- 
independent Lagrangian C = C(F,G), written in terms 
of the two invariants F = F^F^ and G = F» U F* V , the 
associated energy-momentum tensor, defined by 



reduces to 



-AC F F^ a F av + (GjC G -C)g^. 



(4) 



(5) 



p = -C - 4E 2 C F , p = C + -{E 2 -2H 2 )C F . (12) 

o 



III. MAGNETIC UNIVERSE 

A particularly interesting case of the scenario outlined 
in the previous section occurs when only the average of 
the magnetic part is different from zero. We shall call 
such a case a magnetic universe. This situation turns 
out to be relevant in cosmology, because the electric field 
is screened by the charged primordial plasma, while the 
magnetic field lines are frozen ^jjj. In spite of this fact, 
we shall devote some attention to the mathematically 
interesting case in which E 2 = a 2 H 2 ^ in Appendix 1. 

Before discussing the magnetic universe, we shall show 
that when the dynamics of the EM field is given by the 
following series 



C = Y,CkF k 



(13) 



(where k takes values on the integers) the fluid can be 
interpreted as composed of k non-interacting fluids. In 
order to show that this property is indeed valid, we shall 



3 



work with the standard form of the FLRW geometry in 
Gaussian coordinates 



ds 2 = dt 2 



i(t) 2 (dr 2 + X 2 dn 2 ) 



(14) 



The conservation of the energy-momentum tensor pro- 
jected in the direction of the co-moving velocity u M = 8q 
yields 



p+{p + p)6 = Q, 



(15) 



where 9 = uf„ = 3 -. 



Using the values of energy density and pressure given 
in Ean. i|12|) along with Ean. (|13|l we obtain that 



P : 



k 



P 



(16) 



where, for the magnetic universe, we have 
Pk = -c k 2 k H 2k , 



Pk 



c k 2 K H ZK I 1 - ^ 



From the conservation equation (|15(l we obtain 



(H 2 )-+AH 2 



0. 



(17) 



(18) 



The important result that follows from this equation 
is that the dependence on the specific form of the La- 
grangian appears as a multiplicative factor. This prop- 
erty allows us to obtain the dependence of the field with 
the scale factor independently of the particular form of 
the Lagrangian, since Ean. (|18fl yields 



Thus, in order to have an accelerated expansion, the 
source must satisfy the constraint (p + 3p) < 0. In 
terms of the quantities defined in Ean. p2|l . this con- 
straint translates into 



C F > 



C 
AH 2 ' 



(22) 



Any nonlinear electromagnetic theory that admits a con- 
figuration satisfying this inequality will yield acceler- 
ated expansion. Referring now to the particular La- 
grangian given in Ean. l|13fl . Ean. (|2L)fl shows that if the 
fc-component with k < i dominates the evolution of the 
universe, there will be accelerated expansion (assuming 
that p k > 0). 

The conditions to have a bouncing universe follow from 
Raychaudhuri's equation, which states that the expan- 
sion factor must obey the equation 



3 2 KH 



■3p). 



(23) 



In addition to this equation, the energy content of the 
fluid must satisfy the constraint 



1 

p= - i 
H 3 



3e 

^2 - 



(24) 



which comes from the — component of Einstein equa- 
tion in the FLRW framework. Thus two conditions must 
be fulfilled for the bounce to be a minimum of the scale 

factor (p + 3p)b and pb = % 1301 . 

a b 

Having established the conditions for accelerated ex- 
pansion and for a bounce, we shall show next that the 
magnetic universe with the Lagrangian given in Eqn.JSJ) 
satisfies both conditions. 



/ / - 



(19) 



This property implies that it is possible to associate to 
each power k an independent fluid with energy density 
p k and pressure p k in such a way that the correspond- 
ing equation of state, which is preserved throughout the 
expansion of the universe, is given by 



Pk 



Ak 

T -1)P*. 



(20) 



Before proceeding to examine the consequences of these 
fluids in the expansion of the universe, let us turn to the 
analysis of the conditions needed for a bounce and for an 
accelerated expansion. 



IV. CONDITIONS FOR BOUNCING AND 
ACCELERATION 

We shall begin by setting the conditions needed to have 
acceleration. From Einstein's equations, the acceleration 
of the FRWL universe is related to its matter content by 



3l = -l( P + 3 P ). 



(21) 



V. A BOUNCING UNIVERSE 

Let us briefly review the scenario proposed in Q , where 
a nonsingular universe was presented, based on the La- 
grangian 



4 



aF 2 , 



(25) 



where a has the dimension of (length) 4 . Notice 
that Maxwell's electrodynamics can be obtained from 
Ecin. (|25|l by setting a = 0, or recovered in the limit 
of small fields. Using the average process introduced in 
Sect.©, it follows that 



p = ±H 2 (l~8aH 2 ), 
p = -H 2 (l-40aH 2 ). 



(26) 



(27) 



It was shown in Sect. (|III|l that H = H /a 2 , where H is a 
constant. Inserting this result into Friedmann's equation 
through Ean. (|26() leads to 



El 

6 a 2 



1 - 



8am 



(28) 



4 



Since the right-hand side of this equation cannot be neg- 
ative it follows that, regardless of the value of e, for a > 
the scale-factor a(t) cannot be arbitrarily small. The so- 
lution of Eqn. (|28[l is implicitly given as 



t = ± 



a(t) 



dz 



6? 



8aH* 



(29) 



where a(0) = oq. For the Euclidean section, Ean. l29() 
can be easily solved to yield 



a' = Hn\ -t 



8 at. 



(30) 



From Ean. H19(l . the average strength of the magnetic field 
H evolves with time as 



H z = - 



2 i 2 + 12 a 



(31) 



Expression l|30|) is singular for a < 0, since there exist a 
time t = v — 12a for which a(t) is arbitrarily small. For 
a > the radius of the universe attains a minimum value 
at t = 0, given by 



H V8a. 



(32) 



Therefore, the actual value of a m i n depends on Hq, which 
turns out to be the sole free parameter of this model. The 
energy density p given by Ean. l|2t)[) reaches its maximum 
value Pmax = l/64a at the instant t = to, where 



in = vTsTi 



(33) 



For smaller values oft, the energy density decreases, van- 
ishing at t = 0, while the pressure becomes negative. 
Only for values of t such that t <, \/4a, the nonlinear 
effects are relevant for the solution. Notice that the so- 
lution given in Ecin. (|3U|) reduces to that corresponding 
to Maxwell's Lagrangian, namely the scale factor for the 
radiation era, for large t. 

Next we shall show that the Lagrangian given in 
Ean. (|25|l yields a cyclic universe. The qualitative proper- 
ties of the evolution of the scale factor can be obtained by 
interpreting Ean. (|28|l as describing the the motion of a 
particle with constant energy — e and position a(t) under 
the influence of the potential V given by 



V = — 



HZ 
6a 2 



1 - 



8am 



(34) 



It follows from the analysis of V that the scale factor a(t) 
bounces between its minimum <Zf,(t) and its maximum 
a b (t) which are provided by the two real solutions of the 
cubic equation 



ear 



0, 



(35) 



where x 



VI. AN ACCELERATED UNIVERSE 

As discussed in the previous section, a positive power 
of the invariant F in the Lagrangian yields a nonsingu- 
lar FLRW cosmology. In fact, were higher values of k 
present, they would dominate the early evolution. Let 
us see what happens when there are negative powers of 
F in the Lagrangian by analyzing the example given in 
Eqn. @ 9] . We restrict our analysis to the magnetic uni- 
verse and set E 2 = 0, with a residual magnetic field, for 
the reasons already explained. In this case, 



IP 
2 



1 



2 BR 2 ' 



where we have set 7 = — /1 8 [14|. Since H = Hq/o 2 , 
evolution of the density with the scale factor is 



the 



m 1 



2H 2 



(36) 



For small a it is the ordinary radiation term that domi- 
nates. The 1/F term takes over only after a = \J Fl^j /i, 
and grows without bound afterwards. In fact, the curva- 
ture scalar is 



R = T»=p-3p: 



v 

H 2 



Using the density given in Ean. (|36Jl in Ean. l|21|l gives 



3- 

a 



2 o 4 



2H 2 



a 4 = 0. 



To get a regime of accelerated expansion, we must have 



H 2 
a 



^-3f,a 4 <0, 



which implies that the universe will accelerate for a > a Cl 
with 



4 \ 1/8 




3^ 



The examples of positive and negative powers of F ana- 
lyzed in the two previous sections suggest that it is worth 
analyzing a complete scenario constructed by combining 
these cases in a single phenomenological scenario. This 
will be done in the next section. 



VII. A COMPLETE SCENARIO 

We shall analyze here the evolution of the toy model 
generated by the Lagrangian 



(37) 
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FIG. 1: The Lagrangian given in as a function of F. 



with the dependence of the magnetic field on the scale 
factor given by H = Ho/a 2 . In Fig.QJ we plot the La- 
grangian as a function of the value of the field. The 
derivative dL/dF has three zeros in points a, 6, a In these 
points p + p vanishes. In the case of pure magnetic uni- 
verse the value of F is always positive. 

As shown in Appendix III, the theory described by 
Ean. lpfTjl applied to the case of a point charge gives re- 
sults that arc in agreement with observation. In the cos- 
mological case, it implies the existence of three distinct 
epochs according to the value of the field. In the bounc- 
ing era the value of the curvature scalar is small and 
the volume of the universe attains its minimum. The 
energy density and the pressure are dominated by the 
terms coming from the quadratic lagrangian F 2 and are 
approximately given by 



P 
P 



^(l-8a 2 ff 2 ), 



H 2 
G 



(1 -40a 2 H 2 ). 



In the radiative era, the standard Maxwellian term dom- 
inates. Due to the dependence on a~ 2 of the field, this 
phase is defined by H 2 >> H 4 yielding the approxima- 
tion 



H 2 



P ' 



P : 



H 2 
6 



Finally the accelerating era takes place when the uni- 
verse becomes large, and the 1/F term dominates. In 
such a case, 



2 H 2 ' 



6 H 2 ' 



We can analyze the whole evolution using the effective 



potential, given by 

V(a) 



A 



B 



Ca° 



a" a" 

The constants in V(a) are given by 



(38) 



A = AaH, 



o i 



D 



C 



j£_ 

2Hr 



and are all positive. 

The analysis of V(a) and its derivatives implies solving 
polynomial equations in a, which can be reduced to cubic 
equations through the substitution z = a 4 . The existence 
and features of the roots of such equations are discussed 
in [2(| . A key point to the analysis is the sign of D, 
defined as follows. For a general cubic equation 

x 3 + px = q, 



the discriminant D is given by 



D 



We will denote by Dy the discriminant corresponding 
to the potential and Dyi that of the derivative of V. 
From the behaviour of the potential and its derivatives 
for a — * and a — > oo we see that only one or three 
zeros of the potential are allowed. In the case of one 
zero (described by Dy > 0, Dy> > 0), there must be an 
inflection point, and the qualitative plot of the potential 
is given in Fig. [2 The position of the inflection point 




FIG. 2: Plot of the effective potential for D v > 0, D v , > 0. 



depends on the values of the parameters of the model. We 
see that this model is nonsingular for any value of e, but it 
does not display a transition from an accelerated regime 
to a non-accelerated one, so we shall move on to the other 
possibility (given by Dy > 0, Dy> — 0), displayed in 



6 




effective metric given by 



FIG. 3: Plot of the effective potential for D v > 0, D v , = 0. 



Fig. [3J This figure shows the qualitative behavior of the 
potential for typical values of the parameters. Again, the 
model is nonsingular for any value of e. The maximum 
and the minimum of V are always well above the line 
e = 1 so there is no transition in the acceleration for this 
value of e. For e = or 1, the model displays a "coasting 
period" about the transition point. 



VIII. CAUSALITY AND NLED 

Before going to the concluding remarks, we would like 
to call the attention of the reader to a salient feature 
of nonlinear electrodynamics in a cosmological setting. 
Most of our description of the universe is based on the 
behavior of light in a gravitational field. In a FLRW sce- 
nario the existence of a horizon inhibits the exchange of 
information between arbitrary parts of the universe. The 
observation of the high degree of isotropy of the CMBR 
generated a conflict with causality: different parts of the 
universe in the standard FLRW geometry could not have 
enough time to homogeneize. The inflationary scenario 
was engineered precisely to solve this problem. But it 
is worth to mention that nonlinear theories of electrody- 
namics also introduce a new look into causality, which 
we will overview very briefly next. 

Let us take as an example the action [3lj 



S = 



-7 L(F) d 4 x 



(39) 



Lp t y u,v — 4 Lff F^ a F a u . 



(40) 



Notice that this effective metric is generated solely by 
the self-interaction of the electromagnetic field since the 
metric coincides with the background geometry when the 
theory is linear. 

For Lagrangians that depend also of F* , an analogous 
effective geometry appears, although there may be some 
special cases in which the propagation coincides with that 
in dictated by the background metric @ . Another feature 
of the more general case L — L(F,F*) is that birefrin- 
gence is present. That is, the two polarization states of 
the photon propagate in a different way. In some special 
cases, there is also bimetricity (one effective metric for 
each state). In some special cases of nonlinear theories 
(such as Born-Infeld electrodynamics) a single metric can 
appear |g. 

Let us apply the effective metric to the case in which 
the electromagnetic field rests on its fundamental state. 
In the case of Lagrangian 



C 



1 u H 
--F- — 
4 F 



(41) 



the fundamental state is the particular solution 
F 2 = 4/i 8 , 

which corresponds to an energy-momentum tensor equiv- 
alent to a fluid distribution characterized by the condi- 
tion p + p = and generates a de Sitter geometry for 
the background metric 7^„ as seen by all forms of matter 
and energy content, as far as we neglect the gravitational 
influence of such remaining matter and energy. 

For the Lagrangian given above, the effective metric 
tensor takes the form: 



'(>//) 



F 2 ' 



8^ 
F 3 



which in the case of the fundamental state reduces to 



7 (e//) 



= ±J-F" a F Q 



This is a very peculiar and interesting situation: the fun- 
damental state of the theory described by the inverse 
symmetric Lagrangian given in Ean. (|41|l generates a de 
Sitter universe as seen by all existing matter with one 
exception, namely the photons, which follows geodesies 
in the above anisotropic geometry g^.^. 



(7 is the determinant of the background metric). Using 
either the traditional perturbation method or the method 
of the discontinuities devised by Hadamard (see @ and 
references therein), it can be shown that nonlinear pho- 
tons do not move on the light cone of the background 
metric. Instead they move along the null surfaces of an 



IX. CONCLUDING REMARKS 

The standard cosmological model furnishes a rather 
complete picture of our universe, but it does not solve 
yet some important problems. Among these, the most 
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notable are the initial singularity, and the accelerated ex- 
pansion. We have shown here that a specific realization of 
nonlinear electrodynamics yields a unified description of 
the evolution of the universe, in the sense that with this 
single matter source we obtain a qualitative description 
of three of the main phases: bounce (instead of singu- 
larity), radiation, and acceleration. The description of 
the matter content with a single field has the advantage, 
when compared to multi-field models, that in principle 
less parameters and initial conditions are needed. The 
above mentioned features are shown in the plot of the 
scale factor in terms of the cosmological time (see FigQJ. 



such that 



E 2 = a 2 H 2 , 



where < a 2 < 1. Let us assume that the Lagrangian 
is given by the power series in Ean. H13|) . and study the 
relation between the energy density and the pressure. It 
follows from Eans. (jTT|) and (JT3J that 



2a 2 



1-(T 2 

2 a 2 -2 



(42) 




FIG. 4: Qualitative plot of the scale factor as a function of t. 



Notice that in spite of the fact that the fluid as a whole 
would have a time-dependent equation of state, only one 
component of the fluid dominates in each of the epochs, 
leading to an effective equation of state with constant 
parameter. 

We have presented here a toy model which shows that 
nonlinear electrodynamics could be used to build a com- 
plete scenario of the evolution of the universe. We leave 
the matter of the agreement between the values of the 
equation of state parameters used here and those that 
follow from observation to a future publication. Even- 
tual modifications of the model by the addition of new 
terms in the nonlinear Lagrangian, following the same 
one-field phcnomenological approach used here, will also 
be considered, if suggested by observation. 



Setting pk — Afc pu we obtain 

1 2k(a 2 - 2) - 3(ct 2 - 1) 



A fc = 



(J 2 - 1 - 2<7 2 k 



(43) 



Thus, the equation of state depends on both the power 
k and on the partition cr 2 , the dependence on these pa- 
rameters being caused by the nonlinearity of the theory. 
Note that in the linear case (that is, for k — 1) the value 
of the ratio p/p reduces to the c-independent value g, as 
is clear by inspection of Ean. (|43|l . On the other hand, 
when there is "equipartition" , such that a 2 — 1 then, 
independently of the value of the power k it follows that 
\=\- Note that for k < the equipartition is not pos- 
sible, the value of the partition is restricted to the open 
domain < a 2 < 1. 

In two special cases (i.e. k = 1 or a 2 = 0) the depen- 
dence of the field with a(t) is the same. In the general 
case a straightforward calculation reduces the conserva- 
tion equation, Ean. H15(l in the case in which a is a con- 
stant, to the form: 



2fccr 2 



F 1 + cr 2 a 
F + 1 - cr 2 a 



= 0. (44) 



It follows that different parts of the fluid interact, so it 
is not possible to treat the fluid as composed by distinct 
non interacting parts. The solution of Ean. (|44|l provides 
the implicit dependence of the field in terms of the scalar 
of curvature and the constant a 2 . Let us see this by an 
example. In the case of 



F 7 
4 + F' 

the solution of Ean. H44|) is given by 



(45) 



Appendix I: Almost Magnetic Universe (AMU) 

Although the main goal of this paper is the study of the 
pure magnetic universe, it is worth to make some com- 
ments on the case in which there is an electric component 



JP-^TT (^2 _ a -4 = constant . 



The energy density takes the form 



H 2 



3a 2 



1 



2(o- 2 - !)• 



M 8 H~ 



(46) 



(47) 
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and for the pressure we find 



G 



H 2 



ba 2 -7 
6(<J 2 - l) 1 



(48) 



We see from these expressions that although one could 
separate the fluid in two parts by defining p = pi + p 2 
and p = pi + P2 the relations between these quantities 
depend on the equipartition constant. Indeed, we have 



and 



P2 



5a 2 -7 



Q2- 



(49) 



(50) 



3(1 -3a 2 ) 

Thus, there is an interaction between the fluids. 

Appendix II: Born-Infeld Electrodynamics 

Although the final version proposed by Born and Infold 
of nonlinear electrodynamics involves both the invariants 
F and G, we shall use in this section the simplified form 
L = L(F). In this case, 



P 2 1-VX), 



where 



X = 1 



2(i 2 



F. 



(51) 



(52) 



Following Born and Infold, have added a constant term 
into the Lagrangian in order to eliminate a sort of cosmo- 
logical constant and to set the value of the field to zero 
at infinity. We will come back to this question later on. 
Using expressions (|12fl we obtain 



Vx 



i -Vx- 



p 2 



(53) 



Appendix III: Static and spherically symmetric 
electromagnetic solution and the asymptotic regime 

In this article we have presented a modification of 
Maxwell's electromagnetism that can describe in a uni- 
fied way several phases of the evolution of the universe. 
The least we can ask of this modification is that it agrees 
with conventional electromagnetism when the electric 
field generated by a point charged particle is considered. 
For a general nonlinear Lagrangian L = L(F), Maxwell's 
equation for a point charge reduces to 

r 2 L F E(r) = const. (56) 

In the case of the Lagrangian given in Eqn.@ we get 



E(r) 4a 2 £(r) 2 



AE(rf 



q_ 



(57) 



The polynomial in E that follows from this equation can- 
not be solved exactly, but to study the dependence of E 
with r we can plot from Eqn.JS7|) the function r = r(E) 
(see Fig|SJ . Although the plot displays two branches both 
for positive and negative q, the unphysical branches can 
be discarded without further consequences. It also fol- 




FIG. 5: Plot of r = r(E) for a positive charge. The plot for 
negative q is the mirror image of this plot. 



(3 



It follows that 



v i 1H2 



3,= ^(VA-1 



(54) 



(55) 



is a positive-definite quantity. Hence the Born-Infeld La- 
grangian defined in Ean. (|51|l cannot accelerate the uni- 
verse (see however |T^L 



lows from the plot that E goes to a constant value for 
r — > oo. By taking derivatives of Ean. (|57ll it can be 
shown that the function E{r) has no extrema [32j. Hence, 
the modulus of the electric field decreases monotonically 
with increasing r, from an infinite value at the origin to a 
constant (nonzero but small) value at infinity. Ean. l|57(l 
then shows that Eqc — ±-1/7. This situation is akin to 
that in the theory defined by the action 
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It was shown in Q] that the static and spherically 
symmetric solution of this theory does not approach 
Minkowski asymptotically: it tends instead to (anti)-de 
Sitter spacetime. 

Regarding the behaviour of the field for small values 
of r, if we compare the term corresponding to Maxwell's 
case in Ean. l|57|l with the other two, we get that for the 
field to be Maxwell-like we need that 



^7 << E 2 « 



1 



16a 2 



With the explicit dependence for the field given by 
E(r) = q/r 2 , it would be possible to set a value for a 
in agreement with the observation by 



a « 



2 ' 



16q 

where r is a reference value set by the experiment. 

Let us finish this appendix with two comments. For 
very weak fields we can neglect the quadratic term in the 
Lagrangian, and the energy-momentum tensor reduces to 
the form 



In the asymptotic regime we can set E w /i 2 , obtaining 
for the energy-momentum tensor 



2 o 



n 



1 3 



(60) 



which has the form of a cosmological constant. By adding 
an extra term in the Lagrangian (as was done in the case 
of Born- Infold) we could eliminate the residual constant 
field at infinity. Such ambiguity does not arise in the 
case of Maxwell's Electrodynamics due to the linearity 
of the equations. However, for nonlinear electromagnetic 
theories, the geometrical structure at infinity is a matter 
of choice: the field is a constant that can be different 
from zero. Such a property can be translated in a formal 
question: what is the asymptotic regime of the geometry 
of space-time: Minkowski or de Sitter? In linear elec- 
trodynamics the answer to that question is unique, but 
if nonlinearities are present, the possibility of a de Sit- 
ter structure must be considered. In theories in which 
a solution different from zero for the equation Lp = 
exists, such a question has to be investigated combined 
with cosmology [33|. 



T„ 



-Lrj^v — ALp F^a F", 



(58) 



which in the case of a point charge (Foi = E(r)) is 

1 



T ° = Tl 



2E 2 



1 

2FT< 



(E 4 - 3/) , 
(£ 4 + M 8 ) 



(59) 
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